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Hindcasting  of  Wind-Driven  Anomalies  Using  a 
Reduced-Gravity  Global  Ocean  Model 


ABSTRACT 

Global  versions  of  the  Navy  Layered  Ocean 
Model  are  used  to  hindcast  ivind-driven  ocean¬ 
ic  anomalies.  These  versions  are  reduced- 
graxfity  with  the  lowest  layer  infinitely  deep 
and  at  rest,  and  grid  resolutions  of'/°  and 
are  used.  Winds  at  the  1000  millibar  level  from 
the  European  Centre  for  Medium-Range 
Weather  Forecasts  (ECMWF)  are  used  as  forc¬ 
ing  functions  for  the  models  over  the  1981  to 
1989  time  frame.  The  ability  of  the  models  to 
reproduce  wind-forced  anomalies  on  intrasea- 
sonal  to  interannual  time  scales  is  studied  by 
comparing  the  model  solutions  ^vith  various 
observational  data  sets.  These  include  satellite 
altimetry  data,  drifting  buoy  data,  and  island 
and  coastal  sea  level  data.  The  effects  of  vary¬ 
ing  horizontal  and  vertical  resolution  are  also 
detailed.  The  models  are  able  to  hindcast  many 
of  the  wind-driven  anomalies;  the  best  correla- 
ticn  is  found  in  the  tropical  regions  where  the 
oceanic  response  to  atmospheric  forcing  is  the 
most  rapid.  At  mid-latitudes,  mesoscale  flow 
instabilities  are  a  major  source  of  oceanic 
anomalies.  However,  the  resolution  of  the 
global  models  used  here  is  not  adequate  for 
these.  In  addition,  the  reduced-gravity  models 
lack  the  barotropic  mode  and  realistic  bottom 
topography,  which  can  play  an  important  role 
in  the  flow  instabilities. 

INTRODUCTION 

hree  versions  of  a  global  ocean  model  are 
used  to  investigate  the  ability  of  wind-driven 
ocean  models  to  hindcast  and  to  help  identify  a 
particular  class  of  wind-forced  oceanic  anoma¬ 
lies.  This  is  done  by  forcing  the  ocean  models 
using  an  operational  atmospheric  wind  product 
and  then  comparing  the  model  response  to  suit¬ 
able  oceanic  observations  that  are  contempora¬ 
neous  with  the  atmospheric  forcing.  No  assimi¬ 
lation  of  time-dependent  oceanic  data  is  used  in 
the  model.  As  discussed  by  Hurlburt  (1984; 

1987),  the  oceanic  anomalies  of  interest  here  are 
a  direct  integrated  response  to  atmospheric  forc¬ 
ing  on  time  scales  much  longer  than  the  one- 
week  time  scale  of  atmospheric  predictability. 
Therefore,  they  are  insensitive  to  errors  in 
depicting  the  daily  fluctuations  in  the  weather 
but  are  sensitive  to  errors  on  monthly  and  inter¬ 
annual  time  scales.  Other  studies  of  this  class  of 
anomalies  include  Haney  ( 1980)  for  the  mid¬ 
latitude  North  Pacific.  Busalacchi  et  al.  (1983) 
and  Cane  et  al.  (1986)  for  El  Niiio  and  the  equa¬ 
torial  Pacific,  and  Pares-Sierra  and  O'Brien 


(1989)  for  sea  level  along  the  west  coast  of 
North  America. 

The  global  ocean  models  used  here  are 
reduced-gravity  versions  of  the  Navy  Layered 
Ocean  Model,  which  is  a  descendant  of  the  prim¬ 
itive  equation  model  by  Hurlburt  and  Thompson 
(1980)  and  which  has  been  significantly 
enhanced  by  Wallcraft  (1991).  The  model  grid 
covers  the  globe  from  72°S  to  71°N.  Results  from 
model  simulations  using  two  different  horizontal 
resolutions  are  discussed;  they  are  the  y°  model 
(0.5°  lat.  by  0.7°  long.)  and  the  model  (0.25° 
lat.  by  0.35°  long.).  This  is  the  resolution  for 
each  model  variable.  Sensitivities  to  vertical  res¬ 
olution  are  also  presented.  The  models  have 
either  one  or  three  active  upper  layers  while  the 
lowest  layer  is  infinitely  deep  and  at  rest.  These 
are  commonly  referred  to  as  1.5-  and  3.5-layer 
models  because  they  contain  two  and  four  lay¬ 
ers  respectively  but  represent  only  one  and 
three  internal  vertical  modes  for  the  ocean.  The 
mean  layer  depth  is  at  250  m  for  the  1.5-layer 
model  and  135  m,  335  m,  and  550  m  for  the  3.5- 
layer  models.  Large  deviations  from  these  mean 
interface  depths  occur.  The  models  contain  real¬ 
istic  coastline  geometry  and  include  the  deep 
marginal  seas,  such  as  the  Bering  Sea,  the  Sea  of 
Okhotsk,  the  Sea  of  Japan,  the  South  China  Sea, 
the  Sulu  Sea,  and  the  Indonesian  archipelago. 
They  have  horizontal  mixing  in  the  form  of 
Laplacian  friction  and  vertical  exchange  of 
momentum,  mass,  and  heat.  The  models  are 
thermodynamic  and  horizontal  density  gradients 
within  each  layer  can  be  modified  by  advection, 
diffusion,  entrainment,  and  relaxation  to  an 
armual  mean  density  climatology  based  on 
Levitus  (1982). 

Winds  are  the  only  time-dependent 
forcing  for  the  models  and  no  synoptic  ocean 
data  are  assimilated.  They  were  spun  up  from 
rest  to  statistical  equilibrium  using  the 
Hellerman  and  Rosenstein  (1983)  (HR)  monthly 
mean  wind  stress  climatology.  To  drive  the 
model  on  interarmual  time  scales,  1000  millibar 
(mb)  winds  from  the  European  Centre  for 
Medium-Range  Weather  Forecasts  (ECMWF) 
were  used.  The  combination  of  global  coverage 
and  a  relatively  consistent  long  time  series  led 
us  to  choosing  the  ECMWF  winds  over  those 
from  other  operational  atmospheric  centers.  The 
velocity  components  at  twelve-hour  intervals 
were  converted  to  wind  stress  using  a  constant 
drag  coefficient  of  1.5  x  lO-''  and  monthly  aver¬ 
ages  were  formed.  In  forming  the  surface  stresses 


PAPER 


f .  Joseph  Metzger 

Naval  Research  Laboratory 
Stennis  Space  Center. 
Mississippi 

Harley  E.  Hurlburt 

Nava!  Research  Laboratory 
Stenms  Space  Center. 
Mississippi 

John  C.  Kindle 

Naval  Research  Laboratory 
Stennis  Space  Center. 
Mississippi 

ZivSirkes 

Institute  for  Naval 
Oceanography 
Stennis  Space  Center. 
Mississippi 

James  M.  Pringle 

Planning  Systems.  Inc. 
Slidell.  Louisiana 


Cover  Caption 

A  comparison  of  sea  surface 
hei^t  (SSH)  anomalies  in 
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front  the  I(K)()  mb  winds,  no  directional  changes 
were  made.  Howev  er,  global  simulations  forced 
by  the  Et'MWF  wiitds  yield  a  southern  branch  of 
the  South  Equatorial  Curreni  (SEC)  and  a  South 
Equatorial  Countercurrent  (SECC)  that  are  too 
strong  in  comparison  with  observations  and  a 
North  Equatorial  Countercurrent  (NECC)  that  is 
1  to  2'  north  of  its  observed  position.  It  was 
determined  that  a  ntore  realistic  ntean  state 
could  be  obtained  in  the  ocean  model  by  using  a 
hybrid  wind  set.  The  ECMU’F  mean  over  tlie 
HI81  to  1989  time  frame  was  subtracted  from  the 
monthly  means  and  replaced  by  the  HR  annual 
mean.  This  hybrid  wind  set  (EVH)  was  used  to 
extend  the  integration  already  spun  up  to  statis¬ 
tical  equiUbriunt  using  the  HR  monthly  winds. 
Thus,  the  annual  mean  solution  would  still  be 
driven  primarily  by  HR  while  seasonal  and  inter¬ 
annual  forcing  would  come  from  ECMUT.  All 
model  results  discussed  use  the  E/H  wind  forcing. 

Because  the  model  is  basically  non¬ 
eddy  resolving  and  the  model  is  only 
marginally  eddy  resolving,  they  do  not  repro¬ 
duce  realistic  mesoscale  variability  observed  in 
rntyor  current  systems,  such  as  the  Kuroshio 
Extension  or  the  Gulf  Stream.  In  addition,  the 
reduced-gravity  models  lack  the  barotropic 
mode  and  realistic  bottom  topography,  which 
can  play  an  important  role  in  the  flow  instabili¬ 
ties.  These  instabilities  are  the  primary  cause  of 
the  observed  mesoscale  (-50  to  500  km)  vari¬ 
ability.  However,  the  model  is  able  to  reproduce 
directly  wind-forced  anomalies  on  intraseasonal 
to  interannual  time  scales.  The  primary  objec¬ 
tive  of  this  report  is  to  examine  this  capability. 
These  wind-forced  anomalies  include  El  Nino- 
Southem  Oscillation  events,  sea  level  variability, 
the  position  of  fronts  and  currents,  and  curre.it 
transport  variability.  The  oceanic  response  to 
atmospheric  forcing  is  most  rapid  in  the  tropics, 
which  is  where  these  wind-forced  anomalies  are 
predominantly  found. 

The  model's  ability  to  respond  to  large- 
scale  wind-forced  anomalies  is  illustrated  in  the 
figure  on  the  outside  front  cover  of  this  issue 
( Cover  1 );  shown  is  a  global  comparison  of  sea 
surface  height  (SSH)  anomalies  from  the  3.5- 
layer  model  and  data  from  the  II .S.  Navy's 
Geodesy  Satellite-Exact  Repeal  Mission 
(GEOSAT-ERM).  The  GEOSAT-ERM  data  set 
was  obtained  from  the  National  Oceanographic 
Data  Center  and  was  corrected  for  orbit  error 
using  the  Goddard  Elarth  Model  (GEM-T2)  gravd- 
ty  model  (Haines  et  al.,  1990)  and  software 
developed  by  Koblinsky  et  al.  (1990).  The  large- 
scale  orbit  error  was  further  reduced  using  a 
technique  outlined  by  Sirkes  and  Wunsch  (1990) 
All  data  for  this  comparison  have  been  mapped 
on  a  1°  by  1°  grid  using  optimum  interpolation. 

The  upper  panel  on  Cover  1  is  the 
mean  SSH  from  the  model  over  the  first  forty- 


three  repeat  cycles  (Novenvl)er  8,  i;t8ft  to 
November  10,  1988)  of  the  GEOS.AT-ERM.  It  was 
obtained  by  sampling  the  model  along  the  seven¬ 
teen-day  repeat  tracks  of  the  GEOS,\T  ERM  al 
the  same  times  as  the  satellite.  The  large-scali' 
oceiuiic  gyre  .systems  are  clearly  evident  in  al) 
basins,  but  at  subpolar  latitudes  the  model  pro¬ 
duces  a  relatively  flat  SSH  field  within  the  interi¬ 
or  of  the  subpolar  gyres.  This  is  a  result  of  posi¬ 
tive  wind  stress  curl,  which  raises  the  interfaces 
between  the  layers.  To  prevent  the  interfaces 
from  surfacing,  vertical  mixing  is  used.  Most  of 
the  wind  energy  in  these  regions  with  flat  SSH  is 
lost  to  the  infinitely  deep  abyssal  layer  and  the 
active  layers  are  shallow  and  flat.  Thus  the 
model  has  no  significant  geostrophically  bal¬ 
anced  currents  or  anomtilies  in  the  interior  of 
tlie  subpolar  gyres. 

The  middle  panel  of  Cover  1  shows  the 
model  SSH  deviation  from  the  mean  described 
above  for  repeat  cycle  41  (September  20  to 
October  7.  1988),  while  the  bottom  panel  is  the 
SSH  deviation  from  the  GEOSAT-ERM  for  the 
same  repeat  cycle.  The  GEOSAT-ERM  mean  was 
calculated  over  the  same  time  period  as  the 
upper  panel.  Some  striking  similarities  can  be 
seen  between  the  two  anomaly  fields.  These 
include  a  distinctive  narrow  band  of  po.sitive 
anomaly  across  most  of  the  Pacific  Ocean 
between  5°  and  10°N.  Negative  anomalies 
appear  on  either  side  of  it  in  the  central  and 
eastern  basin.  In  the  GEOSAT-ERM  data  one 
might  think  this  band  of  positive  anomaly  is  the 
result  of  water  vapor  error  associated  with  the 
intertropical  convergence  zone,  but  the  anomaly 
also  appears  in  the  model,  which  is  not  subject 
to  that  type  of  error.  Both  tdso  show  relatively 
high  sea  level  in  the  western  equatorial  Pacific 
Ocean  and  the  tilt  of  the  SSH  along  the  equator 
is  similar.  Additionally,  a  positive  anomaly  in 
both  data  sets  is  seen  to  extend  into  the  Indian 
Ocean  off  the  west  coast  of  Australia. 

Similarities  between  th.  model  SSH 
deviations  and  those  from  GEt'SAT-ERM  can  he 
seen  in  other  regions,  but  areas  of  disagreement 
also  exist.  The  best  agre-  .nent  is  at  the  low  lati¬ 
tudes  where  there  is  a  .elatively  rapid  response 
to  changes  in  the  wiiid  fields.  The  worst  agree¬ 
ment  is  in  the  subpolar  latitudes  due  to  the  inter¬ 
face  ventilation  discussed  earlier.  Some  agree¬ 
ment  can  be  found  at  mid-latitudes  even  though 
the  model  lacks  mesoscale  variability.  Note  the 
positive  anomalies  in  the  Kuroshio  current 
region  south  and  east  of  Japan. 

A  similar  comparison  can  be  found  in 
Hurlburt  et  al.  (1992)  for  the  1.5-layer  global 
model.  The  results  for  the  1.5-  and  the  3.5-layer 
models  are  similar,  but  a  significant  difference  is 
found  in  the  eastern  equatorial  Pacific  where  the 
negative  (blue)  SSH  anomaly  along  the  equator 
is  strongest  further  to  the  west  in  the  3..5-layer 
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FIGURE  1.  1981-1989  mean  sea  surface  height  (SSHi  field  overlayed  with  currents  from  the  *  3  5-layer  g'obal  model  Only 
a  Subregion  is  shown,  the  Pacific  Ocean  north  ot  20  S  The  contour  interval  is  5  cm  Shaded  areas  represent  model  land  i  e 
regions  shallower  than  200  meters  (witn  a  tew  exceptions) 


model  as  observed  in  the  GEOSAT-ERM  data. 
This  is  dur-  to  vertical  mixing  in  the  3.5-layer 
model  associated  with  the  equatorial  undercur¬ 
rent  in  the  eastern  Pacific.  With  only  one  active 
layer,  the  1.5-layer  model  cannot  have  an  under¬ 
current 

These  results  show  the  model  is  able  to 
hindcast  observed  wind-forced  oceanic  anoma¬ 
lies,  a  positive  nrite  for  the  atmospheric  forcing 
from  ECMWF  as  well  as  the  ocetui  model.  Given 
accurate  atmospheric  forcing,  the  model  should 
exhibit  intraseasonal  to  interannual  variations 
similar  to  those  seen  in  in  situ  data.  After  a  brief 
di.scussion  of  the  model  climatology,  the  remain¬ 
der  of  the  paper  will  focus  on  the  tropical  lati¬ 
tudes  based  on  the  preceding  results.  Com¬ 
parisons  will  be  made  between  the  model, 
dnfting  buoys,  and  island  and  coastal  sea  level 
stations 


MODEL  CLIMATOLOGY 

The  discu.ssion  of  the  model  climatology  will 
concentrate  on  the  circulation  in  the  Pacific 
( )cean  north  of  20'S,  since  the  remaining  model- 
data  comparisons  are  in  this  basin.  Figure  1 
shows  the  11181  to  Hf8!(  mean  SSH  field  over- 
layed  with  currents  for  the  '/'"  3.5-layer  model. 
With  a  few  exceptions,  the  mean  solution  is  sim¬ 
ilar  for  all  the  models  di.scussed  here. 

The  Itu-ge  cyclonic  subpolar  and  anti- 
cyclomc  subtropical  gyres  au'e  separated  by  the 
Subarctic  Front  at  40  to45"N  The  strongest 
|)ortion  of  the  subtropical  g>Te  is  bounded  on  the 


north  by  the  Kuroshio  Extension  at  32'’  to  3r)°N. 
giving  a  double  frontal  stnicture  separating  the 
two  gyres.  The  southern  half  of  the  subtropical 
gyre  consists  of  broad  westward  return  flow . 
which  intensifies  to  the  west  and  south  as  the 
North  Equatorial  Current  (NEC)  at  10'  to  16‘’N 

Near  the  Philippines  coast  at  approxi¬ 
mately  14’N,  the  NEC  splits  into  the  northward 
flowing  Kuroshio  Current  and  the  southward 
Mindanao  Current  (MC).  This  split  point  varies 
by  about  a  degree  among  the  various  models, 
but  all  are  in  close  agreement  with  the  observed 
dynamic  height  from  the  Navy's  Generalized 
Digital  Environmental  Model  (GDEM )  climatolo¬ 
gy  (Teague  and  Hogan,  1989;  Hurlburt  el  al., 
1989);  GDEM  places  the  split  point  near  14.5''N. 

In  the  model  the  MC  flows  around  the 
Mindanao  Eddy,  centered  at  5°N.  127^E.  and 
merges  with  the  SEC  to  form  the  eastward  flow¬ 
ing  NECC.  The  latter  is  centered  near  5°N  and 
extends  across  the  width  of  the  basin.  The 
model  shows  the  SECC  al  7°  to  lO^S,  but  it  is 
much  weaker  than  the  NECC  and  is  confined  to 
west  of  the  dateline.  The  model  also  accurately 
depicts  the  SSH  tilt  across  the  Pacific  basin  al 
the  equator.  Philander  (1990)  indicates  that  the 
observed  west  to  east  decrease  is  approximately 
50  cm,  the  same  as  the  model. 

A  detailed  description  of  the  model  cli¬ 
matology  in  the  other  ocean  basins  will  not  be 
presented  here,  but  the  model  is  able  to  repre¬ 
sent  the  basic  upper  ocean  features  in  these 
domains.  In  the  Indian  Ocean,  the  reversals  of 
the  Arabian  Sea  and  Bay  of  Bengal  current 


^  .  »  or 

-<141 

ifi-  f 

o.-  ,  r  !• 


!  Dt  3t  C*  But  .  *.'!/ 


4  •  1 


LNSFSCT’ED  3 


MTS  Joiirfinl  •  \  f>l  JG.  N'-  J  • 


systems  that  result  from  the  seasonal  monsoonal 
circulation  are  depicted;  however,  the  complex 
eddy  structure  off  the  coast  of  Sontalia  during 
the  southwest  monsoon  is  not  clearly  seen  in  the 
' model.  These  eddies  are  better  defined  in  the 
M"  model.  The  large-scale  circulation  in  the 
Atlantic  Ocean  is  simulated,  although  the  repre¬ 
sentation  of  the  Gulf  Stream  is  unrealistic.  In 
part,  this  can  be  attributed  to  the  lack  of  the 
thermohaline  circulation  and  associated  deep 
western  boundary  currents  that  play  an  impor¬ 
tant  role  in  this  basin  (Gordon,  1986;  Schmitz 
and  Richardson,  1991).  In  the  tropical  latitudes, 
the  models  depict  the  retroflection  of  the  Brazil 
Current  near  6°  to  8°N. 


UPPER  OCEAN  TRAJECTORY 
HELDS 

he  next  model-data  compari.son  focuses  on 
the  western  equatorial  Pacific  Ocean  where 
the  upper  layer  current  field  from  the  model  is 
compared  to  that  determined  from  triyectones 
of  drifting  buoys  drogued  at  10  m  and  acoustic 
Doppler  current  profiler  ( ADCP)  data.  The 
observational  data  were  obtained  during  the 
Western  Equatorial  Pacific  Ocean  Circulation 
Study  111  (WEPOCS  III  )  and  were  analyzed  by 
Lukas  et  al.  ( 1991);  they  provide  an  excellent 
quasi-synoptic  view  of  the  upper  ocean.  Figure 
2a  shows  the  mean  July  through  September  1988 
velocity  vectors  averaged  on  C  squares.  It 
includes  both  drifting  buoy  and  ADCP  data  and 


FIGURE  2.  (a)  Mean  July-September  1988  velocity  vectors  averaged  over  C  squares;  includes  both  drifting  buoy  and  ADCP  data  (Adapted  from  Lukas  et  al 
1991)  Upper  layer  currents  averaged  over  the  same  time  period  for  (fi)  the  1.5-layer  model,  (c)  the  'f°  3.5-layer  model  and  (d)  the  average  of  two  'U'  3  5- 
layer  model  simulations  Layer  two  currents  averaged  over  the  same  time  period  for  (e)  the  K"  3.5-layer  model  and  (/)  the  average  of  two  /C  3.5-layer  model 
simulations  The  reference  vector  above  (a)  through  (d)  is  1  0  m  s-'  while  it  is  0  5  m  s-'  for  (e)  and  (f) 


has  been  plotted  on  the  model  topography  for 
ease  of  comparison.  The  upper  layer  currents 
from  the  models  averaged  o\  er  the  same  time 
period  are  as  follows;  Figure  2b,  1.5-laycr 
model;  Figure  2c,  3.r)-layer  model;  and  Figure 
2d,  3,f)-layer  model. 

Figures  2d  and  2f  are  an  average  of  two 
3..>layer  model  simulations.  The  second  real¬ 
ization  used  the  end  point  of  the  first  as  its  Initial 
state.  With  increased  horizontal  resolution  and  a 
much  lower  eddy  viscosity  (300  m-  s  '  versus 
1,500  m-  s-'  for  the  3.5-layer  model  and  2,000 
m-  s  ‘  for  the  1.5-layer  model),  flow  instabili¬ 
ties  are  more  pronounced  in  the  '/°  model  com¬ 
pared  to  the  models.  Thus,  given  the  same 
atmospheric  forcing,  a  small  perturbation  in  the 
initial  state  could  cause  the  two  hindcasts  to  dif¬ 
fer  significantly.  In  order  to  lessen  these  flow- 
instability  effects,  an  average  of  the  two 
model  simulations  is  used. 

The  agreement  between  the  WEPOCS 
III  data  and  the  three  model  simulations  is  gen¬ 
erally  good.  All  show  the  southward  flowing  MC 
along  the  east  coast  of  the  Philippines  and  its 
bifurcation  in  the  Sulawesi  (formerly  Celebes) 
Sea,  The  WEPOCS  III  data  indicate  the  west¬ 
ward  branch  is  strongest  in  the  northern 
Sulawesi  Sea  with  some  recirculation  in  the 
south.  The  /i°  model  (Figure  2d)  agrees  well  in 
the  region.  There  is  some  disagreement  between 
the  models  (Figures  2b  and  2c)  and  the 
observed  data  concerning  this  westward  branch. 
The  models  show  westward  flow,  which 
strengthens  to  the  south,  across  the  entire 
Sulawesi  Sea.  Despite  this  disagreement,  all 
three  models  and  the  observed  data  indicate  the 
MC  as  the  primary  water  mass  source  for  Pacific 
to  Indian  Ocean  throughflow  via  the  Makassar 
Strait.  This  agrees  with  the  global  simulations 
described  in  Kindle  et  al.  (1989)  as  well. 

The  remaining  branch  of  the  MC  flows 
southward  a  short  distance,  turns  eastward 
around  the  Mindanao  Eddy  (ME)  (8°N,  129°E  in 
the  observations  and  the  '/°  model)  and  feeds 
the  NECC.  All  the  models  depict  this  with  some 
having  more  success  than  others.  South  of 
Mindanao  the  two  models  (Figures  2b  and  2c) 
have  stronger  southward  flow  than  is  observed 
and  the  center  of  the  ME  is  too  far  south.  The  '/° 
model  (Figure  2d)  agrees  best  with  the  WEPOCS 
III  data  south  of  Mindanao.  This  is  possibly 
because  the  model  is  better  able  to  resolve 
the  complex  geometry  of  this  region.  Note,  both 
Sangihe  and  Talaud  Islands  are  seen  near 
3'^— l^’N.  126°-127'’E.  The  surface  flow  between 
these  two  islands  is  in  agreement  with  observa¬ 
tions  reported  by  Wyrtki  ( 1961). 

The  SEC  crosses  the  equator  along  the 
north  coast  of  New  Guinea  and  is  relatively 
strong  in  both  the  models  and  the  observed  data. 
Near  2°N.  134..5'’E,  the  WEPOCS  III  data  indicate 


a  small  northward  branch  emanating  from  the 
SEC  that  feeds  the  NECC.  Both  models  bifur¬ 
cate  at  this  longitude  while  the  -/i'  model  does 
not.  The  main  branch  of  the  SEC  joins  the  MC 
near  3°N,  128°E  to  help  feed  the  NECC.  North  of 
the  SEC  two  anticyclonic  eddies  appear  in  the 
observed  data.  The  one  centered  near  4°N. 

131°E  is  the  Halmahera  Eddy  (HE).  The  location 
and  shape  of  the  HE  is  better  represented  in  the 
two  V°  model  simulations.  The  model  places 
the  eddy  slightly  too  far  south  and  its  shape  is 
somewhat  elongated.  The  anticyclonic  eddy  to 
the  east  of  the  HE  is  seen  in  all  model  simula¬ 
tions,  but  the  position  and  shape  varies  from  the 
observations. 

The  observed  data  show  the  NECC 
centered  along  5°N.  At  133°E  a  large  meander 
exists  with  a  smaller,  less  distinct  meander  al 
139°E.  All  three  models  agree  well  with  the  lon¬ 
gitudinal  location  of  the  first  meander,  although 
the  amplitude  varies  among  them.  The  ampli¬ 
tude  of  this  meander  is  quite  large  in  the  '/?  1.5- 
layer  model.  With  only  one  active  layer,  all  the 
wind-driven  energy  is  trapped  in  one  layer  and 
stronger  currents  exist.  In  contrast,  the  '/?  3.5- 
layer  model  has  a  smaller  meander.  With  an 
increase  in  the  number  of  layers,  barotropic 
instability  is  decreased  in  this  region  and  the 
energy  is  spread  out  more  among  the  layers; 
weaker  upper  layer  currents  exist. 

Figures  2e  and  2f  show  the  second 
layer  currents  for  the  3.5-layer  and  X “  3.5-layer 
models,  respectively.  These  are  included  to 
show  the  dramatic  change  in  the  current  struc¬ 
ture  between  the  upper  two  layers.  In  the  sec¬ 
ond  layer  the  eastward  branch  of  the  MC  flows 
around  the  HE  and  feeds  a  weaker  NECC.  It  also 
joins  with  the  SEC  along  the  north  coast  of  New 
Guinea  and  feeds  the  eastward  Equatorial 
Undercurrent  (EUC).  Finally,  the  V°  3.5-layer 
model  shows  a  westward  current  near  3°-5°N, 
138°-150‘'E.  which  is  not  present  in  the  upper 
layer.  It  branches  to  feed  both  the  NECC  and  the 
EUC. 

In  an  effort  to  quantify  this  analysis, 
vector  correlations  were  computed  between  the 
WEPOCS  III  data  and  the  model  results  after  the 
model  data  were  linearly  interpolated  to  the  1° 
observational  grid.  On  the  order  of  250  points 
were  used  to  compute  the  correlations.  In  the 
observations,  the  NECC  is  the  dominant,  strong¬ 
ly  time-varying  feature.  Thus,  the  correlations 
are  strongly  affected  by  the  model’s  positioning 
of  this  current.  Those  simulations  that  do  well  in 
positioning  the  NECC  and  the  HE  have  higher 
correlations. 

Table  1  shows  the  vector  correlations 
between  the  July-August-September  (JAS)  1988 
drifting  buoy/ADCP  data  and  the  various  model 
JAS  means  for  all  the  years  of  the  simulations 
(columns  2  to  4).  Here  again  the  y°  3.5-layer 
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TABLE  1.  Vector  correlations  between  July-August-September  (JAS)  1988  drifting  buoy/ADCP  data  from  Lukas  et  at  (1991)  . 

and  the  various  model  JAS  means  (columns  2-4)  Vector  correlations  between  the  JAS  1988  model  mean  from  one 
realization  versus  the  various  JAS  model  means  from  the  second  realization  (column  5)  Model-model  correlations  were 
limited  to  those  poihts  where  observational  data  existed 


Observations  vs  Model  Model  vs  Model 


YEAR 

A  1.S  layer 

'A  3.5  layer 

/•  3.5  layer 

'/•  3.5  layer 

1990 

— 

— 

0.39 

054 

1989 

0  48 

057 

053 

0.60 

1988 

054 

0.54 

062 

082 

1987 

047 

0  48 

0.45 

045 

1986 

0  44 

0  44 

0,47 

043 

1985 

052 

0.53 

0.45 

0  47 

1984 

020 

017 

0.15 

-001 

1983 

017 

0.16 

0.05 

0.08 

1982 

0.47 

046 

026 

0.30 

model  results  are  for  the  average  of  two  realiza¬ 
tions.  Except  for  one  case,  the  JAS  1988 
WEPOCS  III  data  correlate  better  with  the 
model  JAS  1988  mean  data  than  with  the  model 
JAS  means  for  the  other  years.  This  result  also 
applies  to  both  individual  realizations  of  the 
model  (not  shown).  The  '/i°  model  results  for 
1988  are  better  than  either  of  the  '/P  models,  and 
the  'A°  model  discriminates  the  agreement  with 
1988  from  the  other  years  better  than  the  other 
models.  The  one  exception  is  the  3.5-layer 
model,  which  produces  a  slightly  higher  correla¬ 
tion  for  the  1989  JAS  model  mean  than  for  1988: 
in  this  model  the  position  of  the  NECC  in  1989  is 
closer  to  the  observations  than  the  hindcast  for 
1988.  The  large  year-to-year  changes  in  the  cor¬ 
relation  with  the  1988  observations  demonstrate 
significant  interannual  variability  in  the  position 
and  strength  of  the  NECC  and  the  HE.  It  should 
be  noted,  however,  that  part  of  the  interannual 
variability  appears  to  result  from  operational 
changes  in  the  ECMWF  atmospheric  forcing 
products. 

Listed  in  column  five  of  Table  1  are 
vector  correlations  between  the  JAS  1988  model 
mean  from  one  reaJization  of  the  '/f  model  ver¬ 
sus  the  JAS  model  mean  for  the  various  years 
from  the  second  realization.  The  model-to-model 
correlation  calculations  were  limited  to  those 
points  where  observational  data  existed. 
Differences  between  these  two  simulations 
are  attributed  mainly  to  flow  instabilities.  The 
verification  year  has  the  highest  correlation 
(0.82).  However,  the  off-year  correlations  to 
1988  are  very  similar  to  the  off-year  observation 
versus  model  correlations  (columns  2  to  4).  and 
they  exhibit  the  same  degree  of  interannual  vari¬ 
ability.  The  correlations  are  only  slightly- 
improved  in  the  off-years  when  comparing  the 
JAS  mean  from  one  realization  with  the  JAS 
1988  mean  from  the  same  realization  (not 
shown)  Although  the  effect  is  small,  this 


suggests  long-term  anomalies  due  to  flow  insta¬ 
bilities  and/or  long-tenn  memory  due  to  differ¬ 
ences  in  the  initial  state  that  might  influence  this 
region.  The  Pacific  to  Indian  Ocean  throughfiow 
via  Indonesia  is  a  likely  contributor  to  this  long¬ 
term  memory.  Although  the  throughfiow 
exhibits  strong  seasonal  variability,  the  annual 
mean  throughfiow  is  relatively  slow  to  equili¬ 
brate  due  to  the  global  origins  of  this  feature. 
Also  computed  but  not  showTi  were  year-by-year 
correlations  between  the  two  Vf  realizations. 

The  1988  value  of  0.82  shown  in  column  5  has 
the  median  correlation  for  the  years  1982 
through  1990,  which  ranged  from  0.73  to  0.94. 

The  detailed  agreement  among  the 
observations  and  the  three  models  is  somewhat 
surprising,  not  just  because  of  the  complexity  of 
the  region  but  also  because  all  of  the  models 
show  the  western  equatorial  Pacific  as  a  region 
that  exhibits  relatively  strong  flow  instabilities 
as  discussed  above.  These  have  characteristic 
time  scales  of  one  to  four  months.  Therefore, 
three  months  is  a  short  time  to  average  them 
out.  In  addition  to  flow  instabilities  and  short¬ 
comings  in  the  models,  gross  differences  in  sam¬ 
pling  and  errors  in  the  atmospheric  forcing  are 
other  m^or  reasons  for  differences  betw-een  the 
models  and  the  observations. 

ISLAND  AND  COASTAL  SEA 
LEVEL  COMPARISON 

he  Integrated  Global  Ocean  Services  System 
(IGOSS)  Sea  Level  Program  makes  isUmd 
and  coastal  data  available  in  the  form  of  monthly- 
means.  The  time  series  represent  deviations 
from  the  1975  to  1980  mean  .sea  level  computed 
for  each  station;  all  data  have  been  quality  con¬ 
trolled  and  any  large  linear  trends  have  been 
removed.  Sea  level  time  series  from  the  model 
have  been  extracted  at  the  gridpoini.s  closest  to 
the  island  or  coastal  station.  Since  the  model 
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lanti/sea  boundary’  is  the  200  ni  isobath,  the  sam¬ 
pling  points  may  differ  by  up  to  a  few  degrees — 
a  point  to  consider  in  comparing  the  results.  The 
long-term  mean  has  also  been  subtracted  front 
the  model  sea  level  time  series. 

Comparisons  between  the  three  global 
models  and  the  observed  IGOSS  monthly  sea 
level  data  were  conducted  for  twenty-nine  island 
and  coastal  stations  in  the  Pacific  Ocean.  The 
mtOority  of  stations  are  in  the  tropical  latitudes, 
but  a  few  are  located  at  mid-latitudes  (e.g.,  Neah 
Bay,  Washington  and  Crescent  City,  California). 
Figure  3  shows  the  location  of  the  IGOSS  sta¬ 
tions  plotted  on  the  model  domain.  The  compar¬ 
ison  for  Tumaco,  Colombia,  a  coastal  station  in 
the  eastern  Pacific  located  near  2°N,  79°W',  is 
shown  in  Figure  4.  The  upper  panel  (Figure  4a) 
shows  the  IGOSS  data  plotted  with  the  VP 

1.5- layer  model;  the  lower  panel  (Figure  4b) 
shows  the  IGOSS  data  plotted  with  the  VP  and  VP 

3.5- layer  models.  In  this  region  all  three  models 
produce  similar  time  series.  In  general,  the  mod¬ 
els  are  able  to  reproduce  the  interannual  sea 
level  variations  quite  well.  This  is  especially  true 
in  the  example  showm;  the  large  increase  in  sea 
level  associated  with  the  1982-83  El  Nino  is 
clearly  seen  in  all  data  sets.  This  provides  a  mea¬ 
sure  of  confidence  in  the  atmospheric  forcing 
data  and  in  the  ability  of  the  model  to  respond  to 
the  forcing  on  these  time  scales. 

On  seasonal  and  intraseasonal  time 
scales  the  models  have  mixed  success.  At  times 
they  are  able  to  reproduce  both  the  phase  and 
amplitude  of  the  sea  level  signal,  but  frequently 
both  are  in  error.  Generally,  the  variability  about 
the  mean  is  also  weaker  in  the  models  than  in 
the  IGOSS  data.  This  is  probably  partly  due  to 


the  use  of  monthly  ntean  wind  forcing.  Since  the 
Nyquist  frequency  is  two  months,  perhaps  the 
best  to  be  expected  is  for  the  models  to  resolve 
the  semi-annual  time  scale. 

To  quantify  this  comparison,  correla¬ 
tion  coefficients  between  the  IGOSS  and  model 
data  have  been  computed  (Table  2);  the  values 
for  the  ‘4°  model  are  again  a\  erages  for  the  two 
realizations.  If  there  were  tw  o  realizations  from 
the  model  and  the  results  were  also  averaged, 
it  is  expected  that  the  correlations  would  be 
slightly  higher.  However,  the  flow  instability 
effects  would  not  be  nearly  as  pronounced  as 
they  are  in  the  VP  simulations.  The  correlations 
are  generally  good,  however,  there  are  a  few 
cases  where  the  model  does  not  accurately 
determine  the  sea  level  (Noumea  and  Midway). 
In  the  eastern  tropical  Pacific,  the  ‘VP  1.5-layer 
model  represents  observed  sea  level  variability 
with  accuracy  equivalent  to  its  3.5-Iayer  counter¬ 
parts;  exceptions  to  this  are  a  few  stations  along 
the  coast  of  Mexico.  But  the  addition  of  multiple 
layers  (compare  columns  ‘VP  1.5-layer  and  VP  3.5- 
layer)  leads  to  an  overall  improvement  in  the 
statistics.  Most  notable  are  Nauru  and  Tarawa, 
w'hich  are  located  near  the  equator  in  the  Gilbert 
Islands. 

As  showm  in  Figure  2,  a  m^or  equatori¬ 
al  current  system  is  added  to  the  model  when 
the  vertical  resolution  is  increased  from  1.5  to 
3.5  layers.  In  the  second  layer,  the  eastward 
EUC  and  the  westward  currents  that  parallel  it 
to  the  north  and  south  are  added  (Figure  20 
The  3.5-layer  models  also  produce  a  more  realis¬ 
tic  SEC  than  the  1.5-layer  model.  In  the  former, 
the  SEC  extends  a  few  degrees  into  the 
Northern  Hemisphere  while  in  the  latter  the 


FIGURE  3.  Locations  of  the  Integrated  Global  Ocean  Services  System  (IGOSS)  Sea  Level  Program  island  and  coastal  stations 
plotted  on  the  model  domain  _ 
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FIGURE  4.  Monthly  sea  level  data  tor  Tumaco.  ColomOia  (1'50’N.  78=44'W):  a  thirty-dav  running 
mean  has  been  applied  to  tne  model  data  (a)  Observed  Integrated  Global  Ocean  Services  System 
(IGOSS)  data  (heavy  linel  and  the  >:  1,5-layer  model  (thin  line)  (b)  Observed  IGOSS  data  (heavy 
line).  7  3.5-layer  model  (thin  line)  and  the  /)  3  5-layer  model  (dashed  line). 


1981  1982  1983  1984  1985  1986  1987  1988  1989 


broad  NECC  spreads  unrealistically  into  the 
Southern  Hemisphere.  Associated  with  improve¬ 
ments  in  the  3.5-layer  models  is  an  improved  cir¬ 
culation  in  the  Solomon  Sea  and  a  higher  corre¬ 
lation  at  Rabaul. 

An  increase  in  horizontal  resolution 
(compare  columns  3.5-layer  and  '/i°  3.5-layer) 
does  not  significantly  improve  the  correlations. 
Overall,  the  statistics  from  the  V°  3.5-layer  model 
are  equivalent  to  or  slightly  improved  compared 
to  the  3.5-layer  model.  The  higher  resolution 
permits  a  greater  fraction  of  the  variability  to 
occur  in  the  form  of  flow  instabilities  that  are 
not  a  direct  response  to  atmospheric  forcing.  In 
some  cases  this  leads  to  substantial  decreases  in 
correlation:  the  most  serious  degradation  can  be 
seen  at  Honolulu.  This  is  a  result  of  large  anticy- 
clonic  eddies  that  propagate  past  Hawaii  in  the 


retiini  flow  of  the  subtropical  gvre.  1  hese  eddies 
have  unrealistically  huge  sea  level  amplitude  in 
this  model,  but  they  are  not  well  resoh  eii  at 
Nor  are  they  found  at  such  large  amplitude  in 
the  '  (i-layer  Pacific  model  witn  realistic  hot 
tom  topography  discussed  by  Hurlbun  et  al. 
(1992). 

In  order  to  deteniiine  the  model 
response  to  interannual  time  scales,  a  twelve- 
month  running  mecin  was  applied  to  the  model 
time  series  and  to  the  IGOSS  .stations  without 
long  data  gaps,  and  correlation  coefficients  were 
computed  (not  shown).  For  stations  in  the  east¬ 
ern  Pacific,  this  correlation  was  higher  than  in 
lable  2  at  all  locations.  The  average  correlation 
for  the  monthly  time  series  was  O.Gl.  but  the 
average  correlation  for  the  twelve-month  run¬ 
ning  mean  data  was  0.82.  The  highest  correla¬ 
tion.  which  reached  0.91.  was  at  Crescent  City. 
California.  The  stations  in  the  western  and  cen¬ 
tral  equatorial  Pacific  did  not  show  this  marked 
improvement  in  determining  the  interannual 
variability  over  the  1981  to  1989  time  frame.  The 
average  correlations  for  the  tw'elve-month  run¬ 
ning  mean  data  were  nearly  identical  to  the  aver¬ 
age  correlations  for  the  monthly  data.  Analysis 
showed  that  the  later  years  ( 1988  and  1989)  of 
the  time  series  seemed  to  lower  the  interannual 
correlations.  When  recomputed  over  the  1981  to 
1987  time  frame,  the  average  correlations  of  the 
twelve-month  running  mean  data  for  the  western 
and  central  stations  increased  by  approximately 
the  same  margin  as  was  seen  in  the  eastern 
Pacific  stations. 

This  lowering  of  correlations  by  includ¬ 
ing  time  series  after  1987  appears  to  result  from 
an  operational  change  to  the  ECMWT  model- 
data  assimilation  in  1987.  which  increased  the 
amplitude  of  the  winds  in  the  tropics.  These 
results  do  not  imply  a  decrease  in  the  quality  of 
the  ECMWF  winds  due  to  the  operational 
change.  After  an  ac^justment  period,  the 
observed  and  hindcast  SSH  time  series  appear  to 
parallel  each  other  at  least  as  well  as  before  the 
change.  The  decrease  in  correlation  from  includ¬ 
ing  the  period  after  the  change  results  from  the 
relative  length  of  the  time  series  before  and  after 
the  change  and  the  strong  signal  of  the  1982-83 
El  Nino  in  the  time  series. 

CONCLUSIONS 

The  ability  of  the  global  reduced-gravity  mod¬ 
els  to  hindcast  wind-driven  anomalies  has 
been  documented  by  comparing  model  results 
with  various  observational  data  sets.  These 
include  satellite  altimetry,  drifting  buoys,  and 
sea  level  data.  Three  models  using  various  hori¬ 
zontal  and  vertical  resolution  were  studied. 
Because  these  models  omit  the  barotropic  mode 
and  are  either  non-eddy  resolving  or  only 
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TABLE  2.  Correlation  coefficients  between  observed  Integrated  Global  Ocean  Services  System  (IGOSS)  sea  level  data  and  the 
global  models  The  IGOSS  data  are  monthly  averages  and  a  mirty-day  running  mean  has  been  applied  to  the  model  results 


Location 

1/2 

1/2 

1/4 

IGOSS  Station 

(Lat.,  Long.) 

1.5-layer 

3.j-layer 

3.5-layer 

Davao.  Philippines 

(07=05'N.  125'38'El 

068 

066 

068 

Malakal.  Belau 

(U7-'20N.  134=28'E) 

0.40 

046 

043 

Yap  F  S.M 

(09  31'N,  138-08'E) 

036 

0.26 

032 

Truk  Moen,  F.S  M 

(07'27'N.  151“5rE) 

0.61 

0.58 

0.37 

Kwaialem.  Marshall  Is 

(08  44’N,  167  44  E) 

0.54 

0.54 

033 

Tarawa.  Kiribati 

(01"22'N.  172=56'E) 

031 

0.52 

053 

Nauru.  Nauru 

(00’32'S,  166=54'E) 

024 

0.47 

0.46 

Rabaul,  P  N  G 

(04'^12'S,  152°'i1'E) 

075 

083 

0.83 

Penrhyn.  Cook  Island 

(og^ors,  i58=04'W) 

0  71 

0  76 

0.77 

Honiara,  Solomons 

(09=’26'S,  159“57'E) 

0  74 

073 

0.82 

Suva.  Fl|i 

(IS'-OS'S,  178^26'E) 

058 

067 

0.56 

Noumea.  New  Caledonia 

(22'18S,  166=26'E) 

0.25 

0.29 

0.21 

Midway.  USA  Trust 

(28'’13'N,  177'’22W) 

-0  02 

0.01 

0.13 

Honolulu.  HI 

(21°18'N,  157“52'W) 

0.45 

042 

0.03 

Nean  Bay.  WA 

(48°22'N,  124°37’W) 

0.85 

0.86 

086 

Crescent  City.  CA 

(41  45'N.  124°12'W) 

0.67 

071 

075 

San  Francisco,  CA 

(37"48'N.  122  2b  W) 

057 

0.61 

0.61 

San  Diego,  CA 

(32°42’N,  117°10'W) 

040 

0.57 

0.55 

Isia  Cedros.  Mexico 

(28‘'11'N.  115''13'V') 

044 

0.59 

0.62 

Cabo  San  Lucas,  Mexico 

(22°53'N.  109“54W) 

038 

056 

0.56 

Puerto  Valiarta,  Mexico 

(20''37'N.  105°15'W) 

0.32 

0.38 

0.36 

Acapulco,  Mexico 

(16’50'N.  99°55'W) 

0.54 

0.52 

0.54 

Salma  Cruz,  Mexico 

(16’10'N.  95=12'W) 

0.36 

044 

0.44 

Buenave  itura,  Colombia 

(03’54'N.  77°06'W) 

0.56 

056 

0.57 

Tumaco.  Columbia 

(Ot'-SO  N.  78<=44'W) 

0.63 

0,63 

0.66 

La  Libertad,  Ecuador 

(02"12'S.  80°55  W) 

0.64 

0.62 

0.66 

Callao,  Peru 

(12°03'S.  77»09'W) 

066 

070 

070 

Baltra.  Ecuador 

(00"26'S,  gO^I/'W) 

0.62 

0.57 

057 

Santa  Cruz,  Ecuador 

(00°45  S,  90'’19'W) 

0.64 

0.67 

0,66 

marginally  eddy  resolving,  they  lack  much  of  the 
mcsoscale  vanability  associated  with  the  mtyor 
western  boundary  currents.  Due  to  this  limita¬ 
tion,  most  comparisons  were  for  wind-driven 
anomalies  m  the  tropical  latitudes. 

The  generally  good  comparison  of  SSH 
anomalies  between  the  model  and  the  GEOSAT- 
F'RM  indicates  that  the  model  is  accurately 
responding  to  the  large-scale  wind  driving  in  the 
tropics,  a  result  that  owes  credit  to  the  ECMWF 
l(KX)  mb  winds  as  well.  At  the  same  time  it  adds 
credibility  to  the  large-scale  orbit  error  correc¬ 
tion  techniques  used  to  process  the  altimetric 
data.  On  a  finer  scale,  the  models  were  able  to 
lepresent  the  complex  current  structure  of  the 
western  equatorial  Pacific  Ocean  as  indicated  by 
comparisons  with  drifting  buoys  and  ADCP  data 
Finally,  the  model's  interannual  variability  of  sea 
level  at  various  tropical  and  mid-latitude  stations 
compares  quite  well  with  the  observed  data 

The  greatest  improvement  in  the 
model-data  comparisons  came  when  the  vertical 
resolution  was  increased  from  1.5  to  3.5  layers. 
This  a  lows  for  the  formation  of  the  EUC  and  a 
much  more  realistic  SEC.  Increasing  the  horizon¬ 
tal  resolution  from  to  did  not  significantly 


improve  most  of  the  model-data  comparisons 
presented  here.  The  success  of  these  compar¬ 
isons  depended  on  a  relatively  direct  ocean 
model  response  to  atmospheric  forcing. 

Although  it  was  shown  that  increasing  the  hori¬ 
zontal  resolution  yielded  some  improvements  in 
repr  senting  currents  in  the  western  equatorial 
Pacific,  increasing  the  resolution  also  increased 
the  fraction  of  the  variability  due  to  flow  insta¬ 
bilities,  which  are  not  a  direct  response  to 
atmospheric  forcing.  This  tends  to  degrade  con 
parisons  aimed  at  the  wind-driven  response. 
However,  it  is  obvious  that  even  finer  resolution 
and  the  barotropic  mode  are  needed  to  accurate¬ 
ly  model  the  mesoscale  variability  due  to  flow 
instabilities  in  the  mid-latitude  and  subpolar 
regions  of  the  world  oceans.  (Hurlburt  et  al., 
1992;  Hogan  et  al.,  1992,  this  issue.) 
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